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Abstract 

We study the inclusive (e, e'N), (e,e'NN), (e, eV), (e, eVJV) reactions in 
nuclei using a Monte Carlo simulation method to treat the multichannel problem 
of the final state. The input consists of reaction probabilities for the different 
steps evaluated using microscopical many body methods. 

We obtain a good agreement with experiment in some channels where there 
is data and make predictions for other channels which are presently under in- 
vestigation in several electron laboratories. The comparison of the theoretical 
results with experiment for several kinematical conditions and diverse channels 
can serve to learn about different physical processes ocurring in the reaction. The 
potential of this theoretical tool to make prospections for possible experiments, 
aiming at pinning down certain reaction probabilities, is also emphasized. 



PACS: 21.60.Ka, 24.10. Cn, 25.30.Fj, 25.30.Rw 



1 Introduction 

Experiments on exclusive (e, e') reactions in nuclei are becoming more and more com- 
mon in present electron laboratories. Experiments concentrate on (e, e'p), (e,e'NN), 
(e, e'ir), (e, e'nN) etc. reactions which, together with the inclusive (e, e') measurements 
should help us understand better the mechanisms of nuclear excitation by electromag- 
netic probes. Proton knockout reactions at low energies have been used as tools to 
learn about nuclear structure, high momentum components in the nucleus, nuclear 
correlations, etc. 0, 0, f|, Q . Two nucleon emission reactions have been encouraged 
with the hope that they can teach us something about nuclear correlations @, |6], [7|]. 



1 



However, meson exchange currents mechanisms and delta excitation are usually com- 
petitive and, together with complications in dealing with final state interactions, make 
the extraction of information on nuclear correlations a difficult task |J . 

Mechanisms for 2N emission in (e, e') reactions have been developed in H, BJ and 



11 . The analogous reaction of two nucleon emission induced by real photons 



has received parallel attention [O, O, TA\. The approaches of the Pavia 



and Gent groups [10], [Ll[] are suited to deal with low energy photons, below the pion 
emission threshold which opens new channels for 2N emission [O]. Both groups take 
into account final state interactions of the two nucleons with the nucleus: the Pavia 
group by means of nucleon distorted waves with a complex optical potential, while the 
Gent group uses a real nucleon-nucleus potential. The distortion of the nucleon wave 
function by a complex optical potential removes all events where the nucleons collide 
with other nucleons. The procedure is suited to study NN emission leading to specific 
final nuclear states. However, if one is interested in "inclusive" 2N emission, meaning 
that one wishes to compare with an experiment where 2N are detected but one is 
implicitly summing over all final nuclear states, including more than 2N emission, one 
has to find a different scheme. The distortion by a real potential does not eliminate 
the events where there are nucleon collisions. On the other hand, it does not account 
either for events coming from one nucleon emission followed by secondary collisions 
of the nucleon which lead to two or more nucleon emission. Such events have been 



taken into account in ref. |15| although they are insufficient to account for the 2N 



spectra. This sort of final state interaction is also adressed in |16| in order to account 
for contributions to (e, e'pp) from (e, e'pn) followed by charge exchange in a secondary 
collision of the neutron. However, only analogue intermediate states are allowed to be 
excited in that approach and in nuclei with T = there is no contribution from this 
channel. None of the above mentioned approaches deals explicitly with pion production 
and hence they miss the most important channel of 2N emission at energies above 
E 7 = 250 MeV which consists of real pion creation followed by the absorption of 
the pion on its way out the nucleus (indirect photon absorption) [|TJ|]. In ref. [|14| a 



thorough work was done for 7 absorption in nuclei, including all possible many step 
processes: a) photon absorption by a pair of nucleons, b) photon absorption by a trio 
of nucleons, c) (7, it) at a point in the nucleus followed by pion and nucleon collisions, 
d) (7, 7r) followed by pion absorption on its way out of the nucleus, etc. The work uses 
probabilities of photon absorption or pion photoproduction as input, calculated with 
a microscopic many body approach [17]. Then in a second step, it uses a Monte Carlo 
simulation approach to account for the possible reactions which can take place in the 
interaction of the particles produced with the nucleus. The approach is rather succesful 
in reproducing the (7,pn) cross sections ]TJ, [19], although it has some difficulties 



in the (7,pp) channel which experimentally shows a very small cross section. Similarly 
the agreement with the (7, n) JH]] and (7, ttN) []TH| channels is also quite good 



The work of ref. flT] has been recently extended to virtual photons in order to 



study the (e, e') inclusive reaction p^j . On top of several subtleties of virtual photons, 
the extension of the microscopic many body work of |j7j requires the study of the 
quasielastic channel, 7*iV — > N, which is absent in real photons (at photon energies 
bigger than 100 MeV, in practical terms). The work of evaluates the inclusive 
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(e, e') cross section and at the same time all the probabilities for the reactions that 
take place in the interaction of 7* with the nucleons in the first step: (7*, n), (7*, N), 
(7*,JVJV) and (j*,NNN). 

The agreement of the results of ref. [E3| with the inclusive (e, e') data is quite 



good, but since the theory is able to relate the strength of the cross section to different 
physical channels, the theory has a higher potential predictive power than just the 
prediction of the strength of the total cross section and it is the purpose of the present 
paper to exploit this potential. 



In the next section we summarize the technical details which are used in in 
order to obtain the probabilities. 

The present work hence aims to get a complete description of all processes which 
can occur in the (e, e') reaction in nuclei in a range of energies between 100 MeV and 
about 600 MeV of the virtual photon. These are typical energies of the Mainz facility. 
In this way we shall evaluate cross sections for (e,e'N), (e,e'NN), (e, eV), (e, e'nN) 
etc. reactions and will compare them with some experiments in order to interprete the 
data and relate them to the relevant physical mechanisms of the reaction. 



2 The evaluation of probabilities in the first step 



In a previous paper |23[ we carried out a thorough evaluation of the different mech- 
anisms entering the (e, e') inclusive reaction. This allowed us to compare with the 
inclusive cross section d 2 a / dVt'dE' (Q f , E' solid angle and energy of the outgoing elec- 
tron) and also with the longitudinal and transverse response functions. The agreement 
with the data in the regions of the quasielastic peak, the dip region and the A peak is 
good. On the other hand there is also agreement with the longitudinal and transverse 



response functions, particularly with the latest analysis of [2|| taking the world data. 



The study of |23| was done in a covariant way and the response functions are 



associated to the virtual photon self-energy. Thus, 

W L = [ d 3 rImU 00 (q,p(f)) 

vre \q\ J 

(1) 

W T = -^r J d 3 rImU xx (q,p(f)) 
Furthermore, the many body techniques used there, as it was the case before in the 



study of the different reaction channels in pion-nucleus 25 and real photon-nucleus 



reaction fL7| , allow one to separate the different contributions to JmlP 1 ' and associate 
them to the different reaction channels in the first step of the reaction. In order to 
make the statement more explicit, let us briefly discuss the origin and the separation of 
these contributions. In ref. we introduced in a systematic way processes which in- 
volve lplh excitation, 2p2h excitation, 3p3h excitation, lplhlir excitation and 2p2hlir 
excitation. Some of the Feynman diagrams which give rise to ImU^ are drawn in fig. 
1. 
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The imaginary part from a Feynamn diagram is evaluated according to Cutkosky 
rules. The rules [[26[] essentially consist of the following: we cut the intermediate 
states with a line at equal times and when all these states are placed on shell in the 
integrations over the internal variables, this contributes to ImU flu . Placing on shell a 
particle means that we take the imaginary part of its propagator. For instance, for a 
pion 



D(q) ^2ie(q°)ImD(q) = 2 % Q(q°) (-tt) S(q 2 



(2) 



and we see technically the meaning of placing the particle on shell. This way of getting 
ImU flu is doubly rewarding. First, because of the economy versus the equivalent and 
more laborious way (used in fact to prove the Cutkosky rules) of using Wick rotations or 
other alternative methods. Second, because placing some intermediate particles on shell 
is telling us which channel of the total cross section we are considering. Cutkosky rules 
are intimately associated to the optical theorem which relates the total cross section 
to the imaginary part of the forward elastic scattering amplitude. Here it is the same: 
one is relating the imaginary part of the self-energy (also forward) to the probability 
of reaction. Then, the different cuts of intermediate states and their corresponding 
contribution to ImU^, provide us with the physical channels of the reaction and the 
corresponding contribution of these channels to the reaction probability, respectively. 




Fig.l Some many body Feynman diagrams whose imaginary parts contribute to the total 
inclusive (e, e'X) cross section. The dotted lines indicate which particles are put on shell 
when evaluating the imaginary part of each diagram. 

If we look at fig. 1 we can see the different sources of imaginary part represented by 
the cuts of the horizontal lines. Hence, in diagram (a) the cut represents the channel of 
lplh excitation or absorption of the photon by one nucleon. The reaction probability 
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associated to this diagram, or equivalently, its contribution to Wl and Wt is calculated 
in ref. [p5| . Here we would like on top of that to have the distribution of momenta 
of the outgoing nucleon. This can be done by not performing the integration over the 
momentum of the nucleon occupied states in [^3|] since p' N for the outgoing nucleon is 
p + q. However, since we generate events probabilistically one by one, we choose the 
alternative procedure of generating a random momentum from the local Fermi sea, and 
this will generate an outgoing momentum p' N = p N + q (p N = p in fig. 1). If it happens 
that \p N \ < kp(r) (the local Fermi momentum) then the event is Pauli blocked, it is 
dismissed and another event is generated. Thus, we have already the configuration of 
the final state after the first step: in this case just one nucleon produced in the point r 
of the nucleus with momentum p' N . As with respect to having a proton and a neutron 
in the final state, this is trivially done, since the contribution of diagram 1(a) in ref. 
was already splitted into a proton and a neutron induced ones (see eq. (71) of ref. 



H). 

In diagram (b) the cut represents Iplhln excitation or, what is the same, the first 



step in (7*, n). Here we have a probability for this to occur, which is evaluated in |23 



and now one has to generate events with a momentum p' N and p v for the outgoing iV 
and 7r. These momenta are generated according to the weight that they have in the 
evaluation of the probability, which actually involves integrations over these momenta. 
The technical way to do it is using the random procedure which will be shown later 
on. 

The diagram (c) has a cut corresponding to 2p2h excitation. Here the channel 
represented by this cut is two nucleon emission in the first step. Again we will have 
to determine the outgoing nucleon momenta through a random procedure but weighed 
by the momenta distribution implicit in the formulae which evaluate the probabilities 
by integration over these momenta. 

In diagram (d) there is a novelty, peculiar to virtual photons. Indeed, the diagram 
has actually two sources of imaginary part, one of them corresponding to the upper 
cut, where we have a 2p2h excitation and another one corresponding to the lower cut 
where we have Iplh excitation. The lower cut has the same structure as (a) and can 
be cast together using the structure of (a) but with the upper vertex renormalized. In 
this case we have 1 particle emission. However, the upper cut would contribute to 2iV 
emission. 

Finally, the diagram (e) contains two cuts, the upper one which corresponds to 
2p2h excitation, or 2N emission, and the lower one, which appears because we have a 
full (non- static) pion propagator, and which corresponds to having lplhlir excitation, 
or equivalently, IN and In emission. 

These diagrams, and the corresponding distributions of pions and nucleons that we 
get from them, provide the configuration of the final particles after the first step of the 
collision. Thus we have a set of nucleons and pions with definite momenta produced 
at the point r of the nucleus. These particles still have to cross the nucleus before they 
are eventually detected, but in their way out the nucleus they can undergo different 
reactions with other nucleons, thus changing their momenta and leading to the ejection 
of more nucleons. These secondary steps will be followed by means of a Monte Carlo 
simulation method and will be discussed in forthcoming sections. 
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3 The Monte Carlo simulation 



3.1 Generating and tracing the pions 

Let us assume that one pion is produced in a nuclear element of volume, d 3 r. Our sim- 
ulation procedure allows us to decide (depending on the corresponding probabilities): 
first, via which process the pion has been produced and its position inside the nucleus. 
Second, the charge and momentum of the produced pion. 




Fig. 2 Lowest order contribution to the pion production channel. The dot stands for the full 
7* TV — > Nir T-matrix, taken from the model of ref. [p3fl. 



As we have shown in a previous paper [23 the pion production channel has different 
sources: (a) the contribution of the diagrams depicted in fig. 2, (b) the quasielastic 
contribution of the A piece (fig. 3) and, finally, (c) the contribution of the diagram 
depicted in fig. 4 which is related to the (7*,7T7t) channel. Obviously, as we discussed 
[231 , the main contribution corresponds to diagrams (a). 



m 




Fig. 3 Contribution of the A-self-energy in the nuclear medium to the pion production chan- 
nel. 

We evaluate the contribution of the diagrams (a) by using the elementary cross 
section for 7*iV — > Ntt and taking into account nuclear medium effects: Fermi motion of 
the nucleons, modification of the A width in the nuclear medium and the corresponding 
renormalization of the different terms. 
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Fig. 4 Contribution to the 2p2hl7r channel obtained from the 7*iV — ► iV7r7r T-matrix (full 



circle) used in ref. [ 23 1 . 



In ref. [23] we calculate da/dQ'dE', for the (e, e') reaction. We shall refer to it 
as the 7* cross section. In particular the contribution of the ( r y*,ir) channel from the 
diagrams (a) will be given by: 



da„ 



d^dE' 



J d 3 r-0^dn* CM e(cos 9 Max - cos 0cm)||[x 



x2 



n 



AP, dQ'dE'dQ* CM + n[ 



da^*n) x 
' dQ'dE'dn*cM 



(3) 



where j d 3 r is extended to the nuclear volume; J d 3 p corresponds to the integration 
over the nucleon momenta; n p and n n are the proton and neutron occupation numbers; 
6 (cos 9 Max — cos 6 cm) is imposed by Pauli blocking (see ref. p2] ) and fl CM is the solid 



angle of the pion in the CM of the 7*iV frame. The index A stand for the charge of the 
pion. The factor |<#|/|<f| is related to the different flux of particles in the nuclear frame 
and in the nucleon frame (\qi \ is the photon lab momentum in the 7*iV system and 
\q |, the momentum in the nuclear system) and it is actually very close to unity [ p2]| . 

Contributions (b) and (c) mentioned before need at least two nucleons and they 
present a p 2 dependence in the limit of small densities. We calculate these contributions 
following the procedure explained in ref. P31 : 



d a b,c 

dn'dE' 



a \k'\ 



d 3 r f ImW\ bc) L^\ P x 



(4) 



<? 4 \k\ (2tt) 2 

where Pa is the probability related to the charge A of the pion (i.e., the isospin depen- 
dence). For (c) these probabilities are: 



P„(.») = 1 ; P ± ^) = 1 T *JL^. 



(5) 



And for (b) , we consider the probability corresponding to a resonant reaction 
(7*, 7r) with a renormalized pion (because this is the main contribution): 



Pn(7T°) 



2 
3 



Z 
3A 



P_(V~ 



N 
3A 



(6) 
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Once the pion is produced, one must decide which process took place ((a),(b) or 
(c)). If we write: 



a = \- d a j m 

dn'dE' ^ dn'dE' y ' 

j 

where the sum over j stands for the (a), (b) and (c) mechanisms given in eqs. (3), (4). 
The choice of the process is made according to these d 2 aj /dfl'dE' . 

The angular distributions (9, 0) for pions in cases (a) and (b) are generated ac- 
cording to the differential cross section da{p/*N — > TtN)/dVL'dE'dVL* CM which appears 
in eq.(3), while for the case of (c) we consider the momenta distributions given by 
three-body phase space. 

As we can see from eqs. (3), (4), the spatial location of each event is easily generated, 
because the cross section is expressed in terms of probabilities by unit of volume. 

Let us show how the selection of f and angle is made for the mechanisms (a) and 
(b) (in case (c) would be done analogously following the phase space distributions). 
We define for each channel of charge: 

g r (cos9) j f(r)= / d(cos9) g r (cos 9) 

i-2-k 

*r,e(0) / g r (cos9)= / d^t r fi{4>) 



F(r) = [ dr'f(r') 
Jo 

/cos 
d(cos9')g r (cos9') 

T(0) = j* dtftrfW) 







We generate three random numbers x,y,z G [0, 1[ and choose one nuclear radius r 
such that F(r) = xF(oo) and two angles 9 and <fi such that G{cos9) = yG(cos9 = 1) 
and T(4>) = zT((p = 2ir). In order to fix the pion kinematical variables, a random 
momentum is chosen from the Fermi sea. 

Once the kinematical variables are fixed, they are used as input in the computer 
simulation code for pion propagation described in refs. [0, which gives us differen- 



tial cross sections for every pion channel. The mentioned simulation takes into account 
the two possible kind of reactions in the way out of the pion: on the one hand, the pion 
can be absorbed, using its energy in ejecting two or more nucleons. On the other hand, 
the pion can "suffer" quasielastic collisions where the pion gives part of its energy to 
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the nucleus. In this latter case the pion can undergo charge exchange in one unit or 
remain with the same charge. Through double quasielastic collisions one can also have 
pion double charge exchange. Details can be seen in ref. ||25|| . 

3.2 Nucleons: excitation mechanisms 

Nucleons can be excited in different ways through induced reactions by a virtual photon: 

(i) Direct 7*-absorption: 

The virtual photon is absorbed and gives its energy to Af nucleons. For instance, the 
mechanism of fig. 1(a) leads to a primary one nucleon emission, while the mechanisms 
of fig. 1(c) leads to a primary two nucleon emission. Details of this can be seen in ref. 
H. 

(ii) (7*, 7r) reactions: 

Above the electropion production threshold the 7*iV — > ttN' reaction can occur. In 
this reaction (with the exception of a small fraction of coherent electropion production) 
one nucleon is excited. The excitation energy of the nucleon N' is small because most 
of the energy is used in creating the pion. Pauli blocking makes the events with more 
than one excited nucleon not too relevant. 

(iii) (ir, 7r') knock-out : 

Once a tt is produced, subsequent quasielastic collisions N(n, n')N f are possible. 

(iv) 7r absorption or indirect 7*-absorption: 

The produced 7r cannot only be multiply scattered but also be absorbed by Af 
nucleons exchanging virtual pions or other mesons. In this case Af must be greater 
than 1 and, in our approach, we consider events up to Af = 3. 

(v) NN- collisions : 

The number of excited nucleons increases by means of their secondary collisions. 
A nucleon can share its energy with another one in such a way that both remain after 
the interaction with energy greater than the Fermi energy. These secondary nucleons 
will appear mainly in the lower part of the spectrum. 

Obviously, in order to explain the nucleon spectra, one must consider all these 
processes simultaneously. 

Finally, in the next subsection we will discuss the mechanism (v) above. 

We discuss here briefly how the nucleons are generated in the first two mechanisms 
quoted above. Those generated in the (tt, it') and n absorption steps are discussed in 
ref. 0. 

• Direct 7*-absorption 

As we commented before, our theoretical approach allows us to determine which 
is the fraction of direct 7*-absorption corresponding to one, two or three nucle- 
ons. Furthermore, it provides us with the probabilities for the different channels 
in each point of the nucleus. The MC selects randomly -according to the distri- 
bution of probability- which process takes place and where does it happen. In 
section 2 we described how the nucleon momenta and charge are generated for 
the liV induced 7* absorption (see fig. 1(a)). For the case of two and three nu- 
cleon induced 7* absorption the momenta of the outgoing particles will be chosen 
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according to phase space distributions. Their charges will be selected according 
to the isospin dependence of the process. For example, the non-resonant terms 
(called background terms) in the absorption of the virtual photon by two nu- 
cleons, correspond mainly to absorption by pairs proton-neutron (pn) because 
most of them imply the exchange of a charge pion between both nucleons. Hence 
they have a weight proportional to p n p p (p p (p n ) is the nuclear density of protons 
(neutrons)). In some of these reactions, a 7r° is exchanged and they lead to the 
absorption by pairs proton-proton (pp) and neutron- neutron (nn), thus having 
weigths and p 2 n respectively as in the case of the resonant term which has a 
contribution of two-body with an isospin dependence 

1 2 4 1 2 

qP P + QppPn + qPu 

and a contribution of three-body which splits into the four possible channels 
according to the relatives probabilities 



18 8 1 

Y^Pp + ^P 2 ppn + -^PpPn + Y^Pn 

Direct (7*, ir) 

Once the virtual photon momentum is fixed, our code selects a random nucleon 
from the Fermi sea and chooses the charge and the momentum of the outgoing 
pion according to the nuclear cross section (see section 3.1). With this, all the 
properties of the excited nucleon are also fixed. 



3.3 Nucleon propagation 

The nucleons in the nucleus move under the influence of a complex optical potential. 
The imaginary part of the potential is related to the probability of nucleon quasielastic 
collisions in the nucleus (and extra pion production at higher energies which we do not 
consider here) . We consider explicitly these collisions since they generate new nucleons 
going outside the nucleus. As with respect to the real part, we use it to determine the 
classical trajectories that the nucleons follow in the nucleus between collisions. 



As done in ref. [|14| we take as the real part of the nucleon-nucleus potential 



k 2 1 /3 ^ 2/3 

F f 2 



It represents the interaction of a single nucleon with the average potential due to 
the rest of the nucleons. This choice of V(r) means that the total nucleon energy is 
the difference between its kinetic energy and the Fermi energy, €f- 

As with respect to collisions in our Monte Carlo simulation we follow each excited 
nucleon by letting it move a short distance d such that Pd « 1 (P represents the 
probability per unit of length for a quasielastic collision ) . The new position (r ') and 
momentum (p f ) are taken from the Hamiltonian equations as 
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r = f + Sf = f+ pd 
p' = p + 5p; 5p=-WMMl f 

which follow from the Hamilton equations or equivalently from energy and angular 
momentum conservation. 

Our code selects randomly, according to the reaction probabilities which will be 
discussed in 3.4, if the nucleon is scattered or not and, in the case of scattering, what 
kind of process takes place. If no collision takes place, we move the nucleon again. 
When the nucleon leaves the nucleus we stop the process and it is counted as a contri- 
bution to the cross section. If a NN scattering is selected instead, we take a random 
nucleon from the Fermi sea and calculate the initial kinematical variables (P M and s, 
full four-momentum of the nucleon-nucleon system in the nuclear frame and invariant 
energy, respectively). Then, a cos9q M is selected, according to the expression given 
in ref. [ I4| . This expression gives us the correct probability given by da NN /dQcM 



(in the same appendix) plus Pauli blocking restrictions. We take also into account 
Fermi motion and renormalization effects in the angular dependence. We take them 
into account by multiplying each event by a weight factor 

e = A(iv 1 ,iv 2 )^ 7v v^ 

where 1/ X(Ni, N 2 ) is the probability by unit of length that one nucleon N x collides 
with another nucleon N 2 . In average, this factor £ is equal to one. 

Our method assumes that the nucleons propagate semiclasically in the nucleus. 
The justification of this hypothesis for real photons is given in ref. ||14|| . In the next 



section we give some detail on the evaluation of the equivalent NN cross section in the 
medium. 



3.4 NN cross sections 



We are using the parameterization of the NN elastic cross section given in the appendix 
of ref. 
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Due to the fact that for particles of low momenta, the Monte Carlo induces 
large errors, we are not considering collisions of nucleons with kinetic energies below 
50 MeV. 

On the other hand, the reaction probability will change due to the nuclear medium 
effects (Fermi motion, Pauli blocking and medium renormalization). Then, according 
to ref. |14||, the expression for the mean free path (A) of the nucleon is given by 



A 



(iVi 



d 3 p 2 



(2vr 



3 ^P2 



A 



-& N i p (s 



{A- Z) &N 
A K ' 



llab\ 



|Pl 



with P, the full four-momentum of the NN system in the nuclear frame; 
(Pi +P2) 2 ; N-l, the outgoing nucleon and N 2 , the nucleon in the medium. 
Furthermore, 



(9) 



P 2 



a 



N!N 2 



da NlN2 

dtlcM— /(g)e(« - \P.p CM \) 



dtt 



CM 



(10) 
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where CM is the NN center of mass frame and 



x 



29(1 - Id) + — e(|x| - 1) 
d 



(11) 



On the other hand 



P°p CM -e FV G 



(12) 



\P\Wcm\ 

where Pcm is the nucleon momentum in the CM frame and ep is the Fermi energy. 
The function f(q) corresponds to 

1 



/(?) 



' / , I ^ ) V,(<,) 



(13) 



where U(q) is the Lindhard function (U(q) = U^ + Ua) and q is the momentum transfer 
in the nuclear frame. Vj is the spin-isospin effective interaction (details of this can be 
seen in ref. [0]). In these expressions, Q(k— \P-Pcm\) takes into account Pauli blocking 
and f(q), the medium renormalization. 



4 Results 

4.1 (e,e'N) reactions. 

We show our results for the cross section of the (e, e'p) process as a function of the 
missing energy (which is defined as E m = u — T p , the difference between the energy 
transfer and the proton kinetic energy. 

In figs. 5-9 we show results for different kinematical configurations. The agree- 
ment with experiment is fair globally, although some punctual discrepancies appear for 
certain pieces of data. The trend of the spectrum is well reproduced, with the cross 
section increasing with the missing energy in the region of energies studied. The region 
of large missing energy comes from pion production, with or without reabsorption, and 
from secondary nucleons originated in the NN collisions. Some of these mechanisms, 
however, are more visible in the two nucleon emission as we show later. 

Further results can be seen in ref.|J7| which show a very good agreement between 



experiment and the present theory. 

In fig. 10 we show a missing energy spectrum for pn emission at u — 400 MeV. 
We see two peaks, one at low energies and the other at energies around 250 MeV. 
In the figure we have also separated the different sources of contribution. The peak 
at low energies corresponds to the original 7* absorption by a pair of nucleons. We 
include there all events corresponding to the case where the pairs leave the nucleus 
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without further interaction as well as those events where any of the nucleons undergoes 
secondary collisions. We also include the contribution from original three nucleon 7* 
absorption, which in this case is relatively small. 

The broad peak at high missing energies in fig. 10 corresponds to original pion 
production. In this case one has the possibility that the pion is reabsorbed, which 
accounts for the largest part of the strength, while there is a fraction, wich increases 
with the missing energy, corresponding to the case where the pion escapes from the 
nucleus. 

In fig. 11 we show the same figure for pp emission. The features are similar, only 
in this latter case the peak at low energies corresponding to primary 7* absorption by 
a pair of nucleons (mostly pp) with or without FSI is about 5 to 6 times smaller than 
in the pn emission case. On the other hand, the broad peak at large missing energies 
corresponding to primary pion production is about 1/3 the size of the corresponding 
pn peak, much as it happened in real 7 absorption. 
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Fig. 5 The C(e, e'p) differential cross section as a function of the missing energy. E e = 460 
MeV, uj = 275 MeV, \q \ = 401 MeV and 9 P = 0°. Experimental data from [g|]. 
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Fig.6 Same as fig.5 with E e = 478 MeV, u = 263 MeV, 
Experimental data from |1C]. 



303 MeV and 6 p = 38°. 
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Fig.7 Same as fig.5 with E e = 478 MeV, u = 263 MeV, \q\ 



303 MeV and 9 P = 113°. 



Experimental data from |1C]. 
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Fig.8 Same as fig. 5 with E e = 509 MeV, u = 310 MeV, \q \ = 343 MeV and 6 p 



-80°. 
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Fig. 10 C(e,e'pn) differential cross section as a function of the missing energy E m = oj — 
T pl - T p2 . E e = 705 MeV, u = 400 MeV and e = 33.6°. We have separated the different 
contributions: 2N+FSI (main histogram at low miss, energies), 3N+FSI and 17T+FSI (main 
histogram at high miss, energies). 
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Fig. 11 C(e,e'pp) differential cross section as a function of the missing energy. E e = 705 
MeV, lo = 400 MeV and e = 33.6°. We have separated the different contributions: 2N+FSI 
(main histogram at low miss, energies), 3N+FSI and l-zr+FSI (main histogram at high miss, 
energies) . 
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In fig. 12 we show results for the (e, eV) differential cross section as a function of the 
pion momentum for two different angles. This corresponds to pions which manage to 
scape the nucleus without absorption, although they are allowed to undergo quasielastic 
collisions. The pion momenta peak around the resonance region where pion production 
has its maximun strengh. 

In fig. 13 we show a picture in which we separate the primary pions produced into 
those which are absorbed and those which manage to escape as a function of r, the 
radius where they are produced. As can be seen, the pions which are produced in the 
interior of the nucleus are those which have more chances to be absorbed. The results 
are similar for other nuclei, except that the fraction of absorbed pions increases with 
A (see ref. f3|). 
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Fig. 13 12 C: separation of primary pions into those which are absorbed (lower histogram) 
and those which manage to escape (upper histogram) as a function of r. E e = 705 MeV and 
6 e = 33.6°. 

In figs. 14,15 we show the same results but now integrated over r, as a function of 
uj. Figure 14 shows that in 12 C about one fourth of the primary pions produced are 
absorbed in their way out of the nucleus and the rest manages to escape. The situation 
is reversed in 208 Pb, as can be seen in fig. 15, where about two thirds of the primary 

pions produced are absorbed in the nucleus and onlu one third escapes. 
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Fig. 14 Same as fig. 13 integrated over r and as a function of uj. E e = 705 MeV and 
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Fig.15 Same as fig. 14 for 208 Pb. E e = 705 MeV and 9 e = 33.6°. 

In fig. 16 we show the cross section for the emission of a proton and a pion in coincidence, 
as a function of T p . We observe that the energies of the nucleons are relatively low, since the 
pion carries in this case most of the energy. At higher energies, the proton energy distribution 
becomes obviously more spread. 
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Fig.16 C(e,e'p) differential cross section plus a pion detected. E e = 478 MeV, w = 263 
MeV, \q\ = 303 MeV and 6 P = 111 . 
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5 Conclusions 



We have carried out a systematic study of different exclusive channels in the (e, e') reaction, 
like (e,e'N), (e,e'NN), (e, eV) and (e,e'Nn). The procedure used to deal with this multi- 
channel problem in the final state of the (e, e') reaction, has been a Monte Carlo simulation 
method in which the probabilities for the different steps have been evaluated microscopically 
before using quantum mechanical many body techniques. 

The procedure allows one to trace back to the primary step and subsequent secondary 
steps the different contributions to a certain final state channel. For instance, in pn emission 
the pn final state can come from direct knockout of a pn pair by the virtual photon. The same 
process with final state interaction of some of the nucleons, from three body knockout, with 
or without final state interaction of the nucleons. Also, from pion production in the primary 
step followed by pion reabsorption of the pion by pairs or trios of particles. This last process 
shows a strong A dependence in which the fraction of pions which are absorbed changes from 
1/4 to 2/3 from 12 C to 208 Pb. This implies a pion absorption fraction proportional to A a , 
with a ~ 1.3, while in pion absorption experiments, with cross sections close to nR 2 , the 
absorption cross section goes as A a with a = 0.66. This means that the indirect way of 
virtual photon absorption discussed above, with pion production plus reabsorption, is more 
sensitive to the intrinsic probability of pion absorption than the pion absorption cross section 
in pion nucleus collisions. In other words, one can learn more about the mechanisms of pion 
absorption studying indirect photon absorption (with real or virtual photons) than from pion 
absorption in pion nucleus experiments. 

We have compared with a limited amount of experimental data here and found good 
agreement with experiment. Some experimental groups are comparing their data with the 
predictions of the present theoretical approach with good agreement so far. 

The value of the present work can be seen form different points of view. From one side, 
the comparison with data in different kinematical regions and different channels can serve to 
assess the strengh of the different mechanisms discussed here. Conversely, one can use the 
present theoretical tool as a prospective device in order to find the kinematical situations 
which make it more efficient the study of certain mechanisms that one choose to investigate. 



A similar tool for real photons [14| has been used efficiently in both directions [18, 19, 2C] 



and it looks clear that the present study for virtual photons can play a similar role in the 
future in order to extract the maximum information available form forthcoming experiments 
at intermediate energy labs. 
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